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ABSTRACT 

We study the extensively monitored Seyfert galaxy NGC 5548 and compare its nu¬ 
clear emission with models of accretion disk with accreting corona. To obtain the intrinsic 
nuclear spectra from the observed spectra we had to estimate and subtract the contribu¬ 
tion from circumnuclear components such as stars, the Balmer continuum and blended 
Fell lines, and the FC2 extended, featureless continuum. The true nuclear spectra were 
compared with a two parameter model of the accreting disk with an accreting corona, 
described by the mass of the central black hole and viscosity. The model that best htted 
the data was for Mbh = 1-4 x IO^Mq and the viscosity parameter a = 0.033. Such a 
low viscosity parameter was necessary to produce the sufficient amount of X-rays. The 
vertical outflow of mass from corona in the form of wind had to be neglected in our model 
in order to ht into high and low states that NGC 5548 underwent. The model also predicts 
the behavior of the overall opt/UV/X continuum of NGC 5548 during the whole hve year 
campaign. 
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1. Introduction 

The Seyfert 1 galaxy NGC 5548 (z=0.0174) is one of a few active galaxies extensively 
as well as systematically monitored over the last ten years as its nucleus is strongly variable, 
bright and the galaxy is relatively nearby. The amount of data collected for NGC 5548 
actually exceeds all the other monitored AGN: NGC 3783 (Reichert et al. 1994, Stripe 

1994) , Mrk 335 (Zheng et al. 1995) or the BL Lac object PKS 2155-304 (Edelson et al. 

1995) . 

The galaxy has been frequently observed since late seventies in the optical band, in 
the UV with lUE, and by a number of X-ray satellites. In 1988 a systematic campaign has 
been started by AGN Watch team. Optical spectra were taken by ground based telescopes 
every ~ 4 days since December 1988 (Peterson et al. 1991, 1992, 1994). UV variability 
was followed with lUE during hrst eight months of the campaign (December 1988 - August 
1989) with the same frequency sampling (Clavel et al. 1991). The X-rays and UV were 
measured simultaneously by lUE and Ginga satellite 11 times (Clavel et al. 1992). In 1993 
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the galaxy was monitored by HST and lUE in the UV band (Korista et al. 1995). Soft 
X-ray variability was monitored by EUVE satellite (Marshall, Fruscione & Carone 1995) 
as well as ROSAT (e.g. Done et al. 1995). 

The collected data have been used for a number of studies of relative delays in various 
energy bands (optical vs. UV, X-rays vs. opt and UV etc.) as well as between the 
continuum and the broad emission lines. These data offer also an unprecedented possibility 
to study the variability of the nuclear continuum thus offering a chance to distinguish 
between the variety of models able to represent the average overall continuum in a typical 
Seyfert galaxy (e.g. Loska & Czerny 1997). 

However, the continuum models can be tested against the data only if we are able to 
determine the observed continuum emission of the nucleus. It means that the observed 
continuum has to be decomposed into true nuclear emission and contributions from cir- 
cumnuclear components. 

This additional continuum emission clearly consists of the (variable) contribution from 
the Balmer continuum and blended Fell lines emitted by the Broad Line Region (BLR) 
(e.g. Wills, Netzer & Wills 1985) and from the stars located at a distance of a few hundreds 
of parsecs from the nucleus (e.g. Yee & Oke 1978). It was suggested, however, that also 
the third component is present, much bluer than the starlight and rather similar to the 
true nuclear continuum but coming from extended region (e.g. Antonucci 1993, Tran 
1995). Since both this component and the true nuclear emission are featureless, the two 
components are called FC2 and FCl, correspondingly. 

In this paper we attempt to make a decomposition of the optical continuum of 
NGC 5548 derived during the monitoring campaign into nuclear emission and the re¬ 
maining three components listed above. Whenever possible, we use two methods at every 
stage, to ensure that the decomposition does not rely strongly on the adopted method. 
Such a decomposition opens the way to subsequent modeling of the variability of the 
nuclear continuum. 

The outline of the paper is the following. In Section 2 we discuss the contribution of 
starlight to the data and we derive its relative normalization from aperture effects as well 
as from the depth of absorption features. In Section 3 we derive the contamination by the 
FC2 component from the relation between the Hp line and the true nuclear continuum at 
5100 A. Finally, the true nuclear emission in the optical/UV band is presented in Section 
4. We compare it with the predictions of the model of an accretion disk with accreting 
corona. The discussion is given in Section 5. 

2. Contamination by Starlight 

The optical data collected during the monitoring campaign are described in a series 
of papers by Peterson et al. (1991,1992,1994). The data are available to the astronomical 
community due to the courtesy of the AGN Watch team. The spectra are calibrated using 
standard software for a corresponding telescope. However, as different telescopes and 
apertures were used, secondary calibration is needed in order to reduce all the spectra to 
the set A with aperture 5.0” x 7.6”. We performed this reduction following the description 
given by Peterson et al. (1994), i.e. we used values of (f) (a point-source correction factor) 
and G (starlight correction factor at 5100 A) appropriate for a given set of data, calibrated 
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to the mean [OIII] (A5007) flux. We transformed the spectra to the rest frame and we also 
corrected the spectra for the effect of interstellar extinction in the Galaxy. The reddening 
is not significant - we adopted the value of E(B-V) equal 0.05, as inferred from the neutral 
hydrogen column (Burstein & Heiles 1984) in the direction of this galaxy, although even 
smaller value 0.035 was recently used by some authors (e.g. Reynolds 1997). We deredden 
the spectra using the standard Seaton extinction curve (Seaton 1979). It results in the 
increase of the flux at 5100 A by a factor 1.174 and at 1350 A by a factor 1.499. We 
adopted these values in further computations. Spectra prepared in such way were used to 
determine the amount of starlight in the standard set A aperture. 

2.1. Stellar Absorption Lines 

The fraction of light due to stars can be estimated from the ratio of the observed 
equivalent width of stellar absorption lines to their average value in the spectrum of a 
template galaxy, which we take to be the spectrum of the nuclear bulge of M31, after 
Coleman et al. (1980), Table 2. We have used this galaxy as a template since it was shown 
by Wamsteker et al. (1990) that the starlight spectrum of NGC 5548 is similar to the 
starlight spectrum of M31. 

For the purpose of starlight estimation the spectrum n58733h from day 8733 (JD 
minus 2,440,000) from set H of Peterson et al. (1994) was used. As spectra from this set 
have the broadest wavelength coverage (together they cover from 3000 A to 10,000 A), 
more absorption features could be identified. It was also important that the spectrum came 
from the epoch at which the galaxy reached its minimum so the contribution of starlight 
to the total flux was more pronounced. At day 8733 the nucleus was rather dim with the 
total flux at 5100 A of 6.48 x ergs j s j crv? j K . The minimum of F\ during the entire 

1988-1994 campaign was equal to 5.70 x ergsjsjerr?jK (day 8816). 

The n58733h spectrum was normalized by dividing the flux by the local value of con¬ 
tinuum (obtained by fitting a spline to the spectrum). A similar procedure was applied 
to the template starlight galaxy. This method removed broad band trends from the con¬ 
tinuum of these spectra and showed the absorption features more clearly. In Figure 1 the 
normalized spectrum n58733h is presented together with the normalized template galaxy 
M31. 

The spectroscopic indicators of stellar population are such prominent absorption lines 
as: the Gall H (A3968) and Gall K (A3934) line and Nal D (A5890,5896) line which are often 
contaminated by interstellar absorption in the host galaxy. Other absorption features such 
as the Gall G (A4304) and the Mgl b triplet (A5167,5173,5184) are usually contaminated 
by the nearby emission lines (the Gall by H.y and the Mgl by NI) so these lines were not 
used for starlight estimation. The lines that were used are: Fel-|-Gal (A6129,6072) and Fel 
(A7406,7458,7511). In Table 1 we show the percentage of starlight estimated from different 
absorption features and defined as the EW of an absorption feature in the spectrum to 
the EW of the same feature in the template starlight galaxy M31. Next we calculate the 
amount of starlight at 5100 A in absolute units taking into account the spectral shape of 
both the starlight and the combined starlight plus nuclear emission. The value estimated 
from the Fel-|-Gal ( 6100 A) feature is 3.4 X 10 ^^ergs js jcm^ jA. and estimated from the 
Eel (~ 7500A) feature is 3.3 x ergsjs jemS jK. 
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The estimate of starlight contribution based on spectral features does not include 
the possible contribution from featureless FC2 component. It is therefore particularly 
interesting to compare it with the determination of starlight by aperture methods which 
may include most of its contribution if this component is extended. 

2.2. Comparison with Aperture Estimates 

Another estimation of starlight can be made by comparing spectra obtained through 
different apertures. By subtracting a small aperture spectrum from a large aperture spec¬ 
trum we can obtain off-nuclear emission from the circumnuclear ring which should be 
entirely due to starlight. We found a few pairs of observations made with distinctively 
different apertures and separated in time not more than four days, so the variations of 
the flux were relatively unimportant. It was also crucial for the different aperture spectra 
to have large wavelength coverage as to obtain the shape of starlight over a broad wave¬ 
length range. Spectra which best htted the above criteria were those from set F (aperture 
3.2” X 6.4”) and H (aperture 4.0" x 10.0”) which together covered 4500 A to 7000 A. An 
example of an off-nuclear spectrum obtained from the subtraction of spectra from these 
sets (n58765h minus n58765f) is shown in Fig. 2, together with the template galaxy M31 
used in the Section 2.1. The ’’starlight” obtained in such a way is much steeper than 
the template starlight and shows a prominent line. The cause of such steepening is 
the fact that in the Peterson et al. (1994) data the spectra at different apertures were 
obtained with different telescopes and most probably the original reduction of the data 
from the small aperture (the 1.6 m Mount Hopkins observations) was much too simplistic 
(D. Dobrzycka, private communication). 

The available data did not allow for the determination of starlight basing on different 
aperture spectra. It also showed that broad band optical spectra are not calibrated well 
enough to use for comparison with the models. In further Sections we constrain ourselves 
to spectra from set A or fluxes at 5100 A for the purpose of modeling. Because data 
from these sets do not cover the wavelength range of the Balmer continuum and the Fell 
emission lines we do not need to worry about this contribution in the spectra in further 
analysis. 

Much better attempt was made by Romanishin et al. (1995), who using a single 
telescope made direct imaging of NGC 5548. The authors obtained the value of 3.4 x 
ergs/s/cm!^/A. for starlight at F';^(5100A) (for aperture 5”x7.6”) which is identical 
with the value of 3.37 ± 0.54 x ergs/s/cm!^/A. found from comparison of lUE and 

5100 A fluxes by Peterson et al. (1991). This method gave the same value of starlight as 
the method based on stellar absorption features (see previous section). 

The method of estimating starlight based on absorption features is more reliable in 
the case of using high quality nuclear and off-nuclear spectra of the same galaxy obtained 
from the same telescope (e.g. Serote-Ross et al. 1996 for NGC 3516). For NGC 5548 
off-nuclear data are not available and simple use of templates from another galaxy may 
lead to highly uncertain results. 

However, since the outcome of our absorption feature method coincide well with the 
aperture based results of Romanishin et al. (1995) we can safely use the value 3.4 x 
ergs/s/cm"^/A. of starlight in further analysis. 
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3. Contribution from Extended Featureless Continuum FC2 


The monitoring of NGC 5548 showed that the continuum of this galaxy as well as the 
Hp flux varies strongly with time. We use this phenomenon to estimate the amount of the 
extended featureless continuum FC2 within the standard set A spectra. 

The broad Hp line is found to vary in response to continuum variations with a lag of 
about 20 days, which varies from year to year. Figure 3 presents the dependence of Hp 
flux (corrected for Galactic reddening) from optical flux at 5100 A(corrected for Galactic 
reddening and starlight contamination). The Hp fluxes were shifted in time by the mean 
delay of the campaign equal 18 days. 

Apart from the short sequence marked with large circles (dates from 7728 to 7766), 
the line responds almost linearly to the underlying continuum. 

The Hp line consists of components originating in the broad line region and in the 
narrow line region. The component that is sensitive to short-time nuclear continuum 
changes is the component that lies nearest to the central engine, i.e. the H^^^. The 
remaining Hp components originate from the Narrow Line Region (NLR) that lies further 
away (the distance to the NLR is of order of a hundred light years) and is not subject to 
such changes. 

We estimated the value of H^^^ by decomposing the line profile into two components: 
a broad line and a narrow line (see Figure 4). We used for that purpose the spectrum 
n58733a (day 8733) as the broad components during that epoch were weak and the relative 
contribution from the narrow components was high, which allowed for relatively accurate 
decomposition. The value of 1.5 x lQ~^^erg/s/cw? for the narrow component was obtained. 

We compared this value with the value expected on statistical grounds. The amount of 
Hp originating from the narrow line region can be estimated based on paper by Osterbrock 
(1989), who showed that the logarithm of the ratio of Hp to [GUI] flux in the NLR for 
Seyfert galaxies is in the -0.3 to 1 range. The value of 0.35 — 11.1 x 10~^^ergsjsjcm‘^ 
is obtained assuming that the [GUI] flux for NGG 5548, as inferred from large-aperture 
observations (Peterson et al. 1991), is 5.58 x lQ~^^ergjsjcm?. Since the observed angular 
size corresponds to 340 pc per arcsec (z=0.0174, d=70 Mpc for Ho—75 km/s/Mpc) the 
size of the region observed in the standard A aperture is 1.7kpcx2.6kpc which covers the 
whole narrow line region. Therefore the value of H^^^ derived from the decomposition 
of the line originates in the NLR and is well within the expected limits. 

We expect that the line should change in agreement with variations of the nuclear 
continuum Fnud as follows: 


Hp{t) = H^^^ + AFr^uclit) (1) 

However, as the observed flux consists of the variable nuclear part and constant contribu¬ 
tion of the FG2 component, the observed relation is the following: 

Hp{t) = A{Fnuci{t) + Ffc2) + B. (2) 

Gomparing the two formulae we can estimate the contribution from FG2 component 

Ffc2 = {H^^^-B)/A. (3) 
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The values of the coefficients A and B are derived directly from the ht of the straight 
line to the data (Fig.3). The value of the coefficient B is 1.95 x lQ~^^ergjsjcm~‘^ and A 
is equal 0.92 in appropriate units. 

The value of equal 1.5 xlO ^^ergslsjcm^ gives us the contribution of the FC2 

component to the continuum at a level of 0.45 ergjsjcm^. It is about 20 % when 

the source is weak and drops down to a few per cent when it is bright. These values are in 
agreement with the estimate made by Tran (1995), who showed that the FC2 continuum in 
Seyfert 1 galaxies account for only about 4 % of the total featureless continuum (dehned 
as the sum of nuclear continuum and the scattered, polarized ’’featureless continuum”) 
while in the Seyfert 2 for 60%-90%. Therefore, even at its low level, NGC 5548 is not 
truly approaching Seyfert 2 galaxy properties, although the contribution from FC2 is not 
completely negligible either. 

Unfortunately, we cannot say anything about the shape of the FC2 component since 
the available off-nucleus spectra are not reliable (see Section 2.2). Therefore, we have to 
ignore its contribution to the optical/UV spectrum in the next Section. 

4. Comparison of the Intrinsic Nuclear Emission with the Accreting Corona Model 

We use the determination of the intrinsic nuclear emission to test a specihc model of 
an active nucleus. That model consists of an accretion disk surrounded by a corona, with 
both an accretion disk and a corona powered directly by accretion (Zycki, Collin-Souffrin 
& Czerny 1995, Witt, Czerny & Zycki 1997). It describes the spectrum of the nuclear 
emission in the entire optical/UV/X-ray band. The structure of the flow is calculated self- 
consistently at every radius, as described by Czerny, Witt & Zycki (1997), and the emerging 
spectrum is integrated over the disk surface so the stationary model is characterized by 
the global parameters: mass of the central black hole, the accretion rate and the viscosity 
parameter a. 

4 . 1 . High and Low State Models 

We start with detailed modeling of the two spectra representing a high luminosity state 
and a low luminosity state. We choose spectra from set A since these are the most reliable 
(see Sect. 2.2). In Fig. 5 we show examples of a low (n58733a) and a high (n57649a) 
flux spectra. The spectra are corrected for starlight, reddening and transformed to rest 
frame. Unfortunately, they cover only a relatively narrow energy range. The lUE data 
are available for high state only. In order to show an example of the spectrum covering 
the entire optical band but to avoid the problem of starlight subtraction in these not well 
calibrated data we plot (see Fig. 6) the difference between the high and low states from 
set H (n58733h subtracted from n57643h). 

The difference between the two states is considerable. For higher states the blue bump 
becomes more prominent, the Hp line broader and the underlying continuum changes slope. 
The data for high and low states are compared with a model of an accretion disk with a hot 
accreting corona. We adopt the value of the viscosity parameter a equal 0.1, as suggested 
by statistical data for quasars (Czerny, Witt & Zycki 1997). For such a low viscosity the 
radial infall of mater is accompanied by strong vertical outflow which was included in the 
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model, as described by Witt et al. (1997). The best fit to the observed data was obtained 
for the following parameters: mass of the central black hole Mbh = IT x 10®Mq, and 
the following accretion rates: rh=0.05 for the high state and rh=0.03367 for the low state 
(Fig.5). The difference between the high and low states resnlting from onr compntations 
is also plotted against the observed difference (Fig. 6). 

In the case of high state the effect of mass loss does not have an important i nfl uence 
on the spectrnm. However, in the case of low state the effect is profonnd so the exter¬ 
nal accretion rate is only slightly lower than in the high state bnt the differences in the 
bolometric Inminosities and spectral shapes between the two states are considerable. 

The model gives a good fit to the observational data. It also predicts the X-ray flnx 
in the 2-10 keV region. The estimated valne of the 2-10 keV flnx obtained from the model 
for the high state is: 2.04 x ergsjs jcm?. 

This estimate can be compared with a prediction based on observational data from 
Clavel et al. (1992). The observed (nndereddened) UV flnx at A = 1350A (logi/ = 15.35) 
at that epoch was eqnal to b.^2xlQ~^^ergs/s/cm?/K. For Fa( 1350A) of the above valne 
the linear formnla gives the X-ray flnx eqnal Q.Q4:XlQ~^^ergs/s/cm?. 

The model therefore is fainter in X-rays than it is expected. Althongh the linear for¬ 
mnla clearly overpredicts the expected valne of X-ray flnx, as there is a strong flattening 
in the dependence of the X-ray brightness from the UV flnx towards higher UV Inminosi¬ 
ties, still the reqnired valne shonld be abont 5xlO“^^er(7s/s/cm^, as snggested by visnal 
inspection of Fig. 4 from Clavel et al. (1992). 

Higher X-ray Inminosities for given UV flnx are expected from the model if the vis¬ 
cosity parameter is smaller. Adopting a = 0.05 we can fit the high state data with the 
parameters: Mbh = 1-4 x IO^Mq, M = 0.037MQ/yr and predict the 2-10 keV flnx eqnal 
?>.QbxlQ~^^erg/s/cm?, closer the expectations based on observed trends. However, in that 
case the reqnired valne of the mass of the central black hole becomes somewhat larger than 
the valnes between 7 x IO^Mq - 10®Mq which resnlted from estimates based on the stndy 
of the Broad Line Region dynamics (Wanders et al. 1995, Done & Krolik 1996; see also 
Loska & Czerny 1997). 

4.2. Predicted Dependencies on Accretion Rate 

In this Section we compare the overall time behavior of the nnclens with predictions 
of the accreting corona model. We fix the valne of M of the mass of the central black hole 
and calcnlate a seqnence of stationary models parametrized by the accretion rate with 
the spectra covering optical/UV/X-ray band. Therefore, we obtain a theoretical relation 
between the flnxes at 5100 A, 1350 A and the integrated 2-10 keV Inminosity, for each 
M. We compare it with the data from the first year of the observational campaign, when 
the lUE data were available and with the data of Clavel et al. (1992). 

The valne of the mass 1.1 x lO^M© snggested by the fit to the optical data in a single 
high state data set does not reprodnce well the overall trend (see Fig. 7). Therefore, we 
also plot a theoretical line corresponding to the central mass of 8 x IO^Mq and the viscosity 
parameter eqnal 0.1, as before. 

We see that the main trend might be well reprodnced assnming a valne of the mass 
of the central black hole abont 9.5 x IO^Mq which is actnally closer to the valne obtained 
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from dynamical considerations. However, smaller mass models lead to even stronger nn - 
derprodnction of the X-ray brightness than higher mass models (see Section 4.1). 

To inspect the problem of nnderprodnction of X-rays we constrnct a theoretical di¬ 
agram of the dependence between the UV flnx and the integrated flnx between 2 and 
10 keV (Fig. 8). We add the observational points from Clavel et al. (1992). Both families 
of models from Fig. 7 nnderprodnce the X-ray emission and the predicted character of the 
dependence is in contrast with the observational points. 

Therefore we modihed onr theoretical model. We neglected the effect of mass loss 
predicted by the model and nsed the same valne of the accretion rate at all radii for a given 
external accretion rate. This led to a qnalitative change of the shape of the dependence 
between UV and X-ray flnx. We adjnsted the mass of the central black hole to represent 
well the relation between the optical and UV flnx (dashed line in Fig. 7). Additionally, 
by taking smaller valne of the viscosity parameter a we conld adjnst the level of the X-ray 
emission to the reqnired level (see dashed line in Fig. 8). 

Onr modihed model of the accreting disk with accreting corona is now consistent with 
the behavior of all three spectral bands: optical, UV and X-ray, remaining still a two- 
parameter model, described by the central mass and the viscosity. In the case of NGC 
5548 these parameters are: Mbh = 1.4 x 10^ Mq and a = 0.033. 

The spectral range in the optical data withont calibration problems (set A) is too 
short to attempt the distinction between the models with and withont mass loss. The 
model withont mass loss, Mbh = 1.4 x 10®Mq, a = 0.033. and accretion rate m = 0.009 
cannot be visnally distingnished from the model given in Fig. 5 (i.e. with mass loss, 
Mbh = 1.1 X lO^M©, a = 0.1. and accretion rate rh ~ 0.034). 

We see that the short seqnence of points marked with large circles (dates from 7728 
to 7766) in Fig. 7 follows approximately the overall trend of the continnnm variations 
althongh the Hp line response was nntypical. Nevertheless, simply on the basis of Fig. 7 
we cannot tell what cansed the departnre of those points from the main trend shown in 
Fig. 3. Two possibilities are open. Either the line flnx was too high (i.e. the BLR did 
not respond to the decrease of the nnclear emission althongh the considered period was 
longer than the nsnal time delay), or the continnnm flnx seen by ns was too low. This 
second possibility wonld happen if some intervening optically thick matter crossed onr line 
of sight. 

Althongh the model of accreting corona is able to reprodnce well the broad band time 
behavior of the nnclens, two problems remain. The hrst is the qnestion, why the mass 
loss predicted by the model is possibly too large. The second is that the favored valne of 
the central mass is somewhat larger than snggested by the dynamical stndies of the Broad 
Line Region, althongh the discrepancy is not stringent. 

5. Discussion and Conclusions 


5.1. Data Requirements 

The character of variability of the Seyfert galaxy NGC 5548 is very complex. It 
means that the observational data not only have to be dense in time, in order to cover well 



the short time scale behavior but also have to extend for long times (years) to catch the 
important long time scale trends. This general problem shows up in a number of particular 
difhculties. 

A few observational points do not always follow the general trend. Two observations 
made when the emission lines in NGC 5548 were exceptionally weak suggested that the 
almost complete disappearance of the lines is a break off from the usual behavior of the 
source (lijima 1992, lijima, Rafanelli & Bianchini 1992, Loska, Czerny & Szczerba 1993). 
Actually, the behavior of the source during that period followed well the general trend, as 
shown by the continuous frequent sampling (Peterson et al. 1994; see also Figure 3). 

On the other hand, there seem to be periods when the source does not follow the 
typical variability pattern. During the period of about six weeks in July/August 1989 
(days 7722 to 7766) the continuum was much fainter than expected on the basis of the 
Hp intensity. Since the line delay (about 18 days) is much shorter than six weeks the 
effect cannot be caused by the lack of answer of the Broad Line Region to the change of 
the continuum. It may suggest some kind of geometrical effect which decreased the level 
of continuum seen by an observer but not by BLR. However, it may also indicate that 
the coherent response of the BLR require more time than indicated by the measured time 
delay. In order to avoid problems, we excluded the data collected during that period from 
our modeling. 

In our comparison of the data with the model the observational point in X-rays ob¬ 
tained in 1985 plays a crucial role since the amplitude of the variation reached exceptionally 
large value at that time. Such events are rare in a sense that large variations happen in 
time scales of years but their inclusion in the study is crucial from theoretical point of view 
so the good data have to cover such long time scales as well as short time scale variability. 

The presence of the broad and continuous range of time scales in the problem is best 
seen in the study of X-ray power spectra which are basically of a power law shape in all 
AGN (McHardy & Gzerny 1987, Lawrence et al. 1987; see also Gzerny & Lehto 1997). It 
is not clear whether the long term variability (in the time scale of years) is actually of the 
same nature that the best studied variability in time scales of hundreds to thousands of 
seconds but we cannot a priori reject such a possibility and the large amplitude variations 
like the one mentioned above should be included in any model unless there are arguments 
against it. 

The currently available data on NGG 5548 are already good enough to allow interesting 
study of the time dependent properties but extension of the observational campaign is 
essential for better understanding of the behavior of accretion flow. 

5.2. Determination of Intrinsic Spectrum of NGC 5548 

The results of Section 2 show that the level of starlight at 5100 A (in the aper¬ 
ture 5.0” X 7.6”) equal to 3.4 x 10“^^ ergs/s/cm^/A is determined reliably. However, 
the spectral shape of the starlight does not seem to be described accurately. Separate 
spectrophotometric observations with a good telescope are necessary in order to solve this 
problem. 

The lack of good starlight determination in the entire optical band limits the testing of 
the model against the data basically to the normalization at 5100 A . It is not a big problem 
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if the entire optical/UV/X-ray range is stndied bnt nevertheless it makes impossible to nse 
the entire optical range for hne tnnning. 

The snbtraction of the FC2 component is another problem which shonld be dealt with 
while hne tnnning the model since the level of the FC2 at 5100 A is np to 20% when the 
sonrce is exceptionally faint. Fortnnately, the sonrce is brighter most of the time so the 
FC2 contribntion conld have been neglected in the present analysis of the applicability of 
the model. 


5.3. Model Advantages and Restrictions 


The variations of the intrinsic continnnm emitted by the nnclens of NGC 5548 in the 
optical/UV/X-ray band were snccessfnlly reprodnced nsing a simple model of an accretion 
disk with an accreting corona by hxing the mass of the central black hole at 1.4 x IO^Mq, 
the viscosity parameter at 0.033 while varying accretion rate from ~ 0.009 to ~ 0.04MQ/yr. 

The model reprodnces well both - the hardening of the spectrnm with an increase 
of the Inminosity and the weak increase in X-ray Inminosity with the increasing UV flnx 
when the sonrce is bright. The disk emission does not extend into soft X-rays so any 
traces of the soft X-ray excess in this sonrce (Walter & Fink 1993) have to be attribnted 
to the atomic processes connected with reprocessing of X-rays by the disk snrface (Czerny 
& Zycki 1994). We neglected those processes in onr model assnming that the disk albedo 
is eqnal to zero. 

The version of the model which corresponds well to the behavior of the sonrce is 
the model in which the vertical ontflow of the mass from the corona in the form of the 
wind is neglected. The theory nnderlying the description of the ontflow will have to be 
reconsidered. Apart from that, the model is a very promising one and shonld be fnrther 
tested and developed. 

The most important topic is the nnderstanding of the natnre of the variations and the 
character of the global response of the disk/corona system. 

In onr present approach we simply nse a set of stationary models characterized by a 
single accretion rate at all disk radii at any given moment. This is certainly an oversim- 
plihcation. The variations in the optical/UV band are noticeable starting from the time 
scales of one - two days (Korista et al. 1995). The main contribntion to the spectrnm 
comes from the inner parts of the disk, at radii r ~ 10 Rschw The thermal time scale of 
that region depends on the viscosity parameter a and is eqnal 


'^th 


0 . 2 ( 


0.033 M 

a ’^IAxIO^Mq’ 


days 


and a more accnrate valne of the two-folding time scale at 1350 A (independent from the 
assnmption abont the disk radins dominating the emission at a given freqnency) is eqnal 
1.0'^ for the mean Inminosity level (Siemiginowska & Czerny 1989). Therefore, if there 
are any changes of the corona strnctnre the thermal strnctnre of the disk can follow these 
changes over the time scales of days. 

The viscons time scale is longer than the thermal time scale by a few orders of mag- 
nitnde (see e.g. Siemiginowska, Czerny & Kostynnin 1996). The disk therefore is never 
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stationary and onr assnmption is not jnstified. On the other hand this assnmption is the 
only reasonable approach since we cannot calcnlate the actnal time behavior of the disk 
becanse of the lack of proper nnderstanding of the disk/corona conpling. The disk/corona 
system described by the model is thermally nnstable (see Fig. 7 of Witt et ah 1997). The 
sndden increase of the corona activity can therefore resnlt even in enhanced disk response, 
as reqnired, bnt this response does not have to satnrate at the valne of the flnx radiated 
by the disk which is necessary for the corona to reach the new temporary thermal balance. 
However, if we do assnme that the thermal disk rnnaway satnrates at the appropriate 
emitted flnx then the disk/corona strnctnre wonld be observationally the same as in the 
case of enhanced accretion rate since the disk/corona conpling does not depend on cold 
disk strnctnre (snrface density etc.). 

The applicability of the stationary model to the description of the time evolntion 
snggests that snch a satnration actnally takes place bnt its mechanism is nnknown. The 
additional disk/corona conpling might be based on magnetic held and the description of the 
disk and corona how by the viscosity coefficient a is too simple to reprodnce well the time 
behavior, inclnding stability, althongh may approximate the qnasi-stationary sitnation. 
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FIGURE CAPTIONS 

Fig. 1. The spectrnm of NGC 5548 (n58733) (thin line) and the template starlight 
galaxy M31 (thick line), both normalized as described in Section 2.1. 

Fig. 2. The off-nnclear spectrnm of NGC 5548 calcnlated by snbtraction of the observa¬ 
tion n58765f (apertnre 3.2” x6.4) from n58765h (apertnre 4.0” xlO.O”). It differs consider¬ 
ably from the overall shape of the starlight template (given with arbitrary normalization). 

Fig. 3. The dependence of the H/^ flnx measnred in nnits of 10“^^ erg/s/cm^ on 
the continnnm flnx at 5100 A measnred in nnits of 10“^® erg/s/cm^/A (after starlight 
snbtraction). Both flnxes were dereddened. The data cover the hrst fonr years of the 
campaign (Peterson et al. 1994). Points marked by large open circles come from the fonr 
week seqnence in Jnly/Angnst 1989 and were not inclnded in the linear ht to the trend 
(dashed line). 

Fig. 4. Decomposition of the Hp prohle into broad and narrow components in the data 
8733, taken when the sonrce was faint. A power law continnnm was adopted, normalised 
ontside the spectral band dominated by strong emission lines. 

Fig. 5. A low (n58733a) and a high (n57649a) states are shown, starlight snbtracted 
and dereddened and the lUE data points corresponding to the high state. The htted 
models are for Mbh = 1-1 x IO^Mq,^ = 0.1 and m=0.05, m=0.03367 for the high and 
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low state spectra, respectively. Outer disk radius was put at 120 Rschw Fluxes are given 
in erg/s/cm^/Hz. 

Fig. 6. A broad band difference between two states (n58733h subtracted from n57643h) 
is shown (in order to avoid stellar subtraction problem), together with the difference be¬ 
tween two theoretical spectra from Fig. 5. 

Fig. 7. The relation between the optical and UV flux during the first year of observa¬ 
tional campaign given by Clavel et al. (1991) (squares) and during the fifth year given by 
Korista et al. (1995) (triangles), after dereddening of the data. Points from July/August 
1989 period with exceptionally high Hp are marked by large circles. Continuous lines show 
the dependence predicted by the model assuming the viscosity parameter equal to 0.1, 
for two masses of the central black hole. The dotted line marks the predictions for mass 
loss excluded, black hole mass equal to 1.4 x 10 ®Mq and viscosity parameter equal to 
0.033. Optical flux is given in units of 10“^^ erg/s/cm^/A, and UV flux in units of 10“^^ 
erg/s/cm^/A. 

Fig. 8. The relation between the UV flux and the integrated 2-10 keV flux. Data 
points are from Clavel et al. (1992), without dereddening. The brightest point represents 
the measurement made in 1985. Lines mark the same sets of models as in Fig. 7, but 
reddened by Galactic extinction. UV flux is given in units of 10“^^ erg/s/cm^/A and 
X-ray luminosity in units of 10“^^ erg/s/cm^/Hz. 


TABLE 1 


starlight feature 

% of starlight 

remarks 

Fel(4923, 5018,5169), Mgl+MgH 

5.8 % 

Mgl contaminated by NI 

Nal(5890,5896) 


emission 

Fel+Cal(6072,6129) 

52.8 % 


CaH(6934) ,FeI(7008,7109,7157) 

12.9 % 

CaH contamination by ISM 

Fel(7459,7511) 

64 % 


Call (8498,8542) 

11 % 

contamination by ISM 
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